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Abstract 

The algorithm of Micali and Vazirani [J] finds a maximum cardinality matching in time 
0{y/nm) if an efficient set-merging algorithm is used. The latter is provided by the incremental- 
tree set-merging algorithm of [5] . Details of this application to matching were omitted from [5] 
and are presented in this note. 

1 Introduction 

Micali and Vazirani implemented the “shortest augmenting path” approach to find a maximum 
cardinality matching in time 0{y/nm) [3]. This time bound assumes a linear-time algorithm for 
set-merging is used. It is tempting to conjecture (as was done in m) that the special structure of 
blossoms achieves this goal automatically. However this claim remains unsubstantiated. Instead 
the incremental-tree set-merging algorithm of [2] can be used to represent blossoms and achieve 
the desired time bound. The details of this application to matching were omitted from [2] and they 
are presented in this note. For accuracy we use the recent paper [5] as the source of all details for 
the cardinality matching algorithm. We assume familiarity with that paper0 

We will present an implementation of the Micali-Vazirani algorithm that achieves time 0{m) 
per phase on a RAM. This easily implies the desired time bound 0{y/nm). The issue is the use of 
set-merging to represent blossoms. Towards this end we divide each phase into two subphases. The 
first subphase assigns levels to vertices. In this subphase all free vertices are essentially equivalent. 
This allows the search structure to be modeled by a tree, and allows use of incremental-tree set 
merging to model blossoms. The second subphase finds augmenting paths. Here we choose bridges 
for DBFS in a depth-first manner. This gives a linear order to newly created blossoms, allowing 
set merging to be handled with a stack. (This stack-based set-merging is used in [T] to compute 
strong components efficiently.) The first subphase requires a RAM, the second subphase runs on 
both RAM and pointer machine. 

The rest of this section gives background from [5]. Then the set-merging algorithm is presented, 
the first “assignment” subphase in Section [2] and the second ’’augmentation” subphase in Section 

El 

The entire discussion uses this terminology: A vertex is a vertex of the given graph G. A 
blossom is a blossom from the previous search level (i.e., tenacity < 2f -|- 1 during search level f); a 
base is the base vertex of such a blossom, or a vertex of finite evenlevel and empty blossom. Thus 
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at the start of a search level, the vertices of G with finite evenlevel are partitioned into blossoms 
and bases. 

As in [21 Definition 29] every vertex v has a value bud*(u). A petat is the set of all vertices 
with finite evenlevel that have the same bud*. Assume a blossom has bud equal to its base, so 
when a search level begins, each blossom plus its base is a petal*. At any point in time the petal*’s 
partition the set of vertices with finite evenlevel. A free petal* is a petal* whose bd* is a free vertex. 
We say bud*(P) = 6 if P is a petal* corresponding to bud* b. 

Note that bud*(u) = u iff u is the bud of a petal*, or (as in M) V is not in any petal*. 

The following lemma is adapted from [5]. We will use parts (i) and (iii)', we include (m) for 
completeness. 

Lemma 1.1 Consider an arbitrary petal* P* with bud b. Let v be an arbitrary vertex in P*. 

(i) Every minlevel(u) path p goes through b. Its subpath from v to b, denoted minlevel(u; b) (like 
the notation of is eontained in P*. 

(ii) The same holds for a maxlevel(u) path p whose (unique) tenacity(u) bridge has been pro- 
eessed by MAX. 

(iii) P* is precisely equal to the union of the paths minlevel(u, b), v & P*. In faet we can restrict 
the union to vertiees v in the support of P*, assuming we add in the petal* whose hud* is v. 

Remarks: Note that (ii) need not hold for a maxlevel path using a bridge that has not been 
processed by MAX (i.e., the petal has not been expanded to its complete blossom). 

Also note that parts {i)-{ii) imply the Micali-Vazirani algorithm handles augmenting paths 
correctly: Suppose P* of the lemma is on an augmenting path; clearly it contains b. Suppose for 
contradiction that v G P* is on another ap (ap = shortest augmenting path) Q that does not con¬ 
tain b. Part (z) shows Q cannot contain a minlevel(u) path. [H Theorem 26] shows a maxlevel(u) 
path M in Q uses a bridge of tenacity < Im- So part {ii) applies to M. Thus M contains b and 
cannot be in Q. 

Proof: Let t be the tenacity in the current search level. 

Define a subgraph H{P*): Start with the graph in which P* is formed (i.e., previous augmenting 
paths, if any, have been deleted). Then remove every bridge that has not been processed by MAX. 
Apply la Lemma 30] to P*. If P* has no tenacity t vertices the lemma shows P* — b is a blossom 
of tenacity t — 2. In the opposite case P* — 6 is a blossom of H{P*) of tenacity t. Since H{P*) 
has IjYi, > t, [a Theorem 7] applies. It shows that as in parts (z) and (zz), each minlevel(u) and 
maxlevel(z;) path M passes through b.. 

[5[ Definition 17] of blossom, plus a simple induction, show M is contained in P*. This com¬ 
pletes the proof of (z) and (zz). It also proves {iii). □ 


2 Assignment Subphase 

Each phase begins in the Assignment Subphase. This subphase assigns levels, according to the 
procedures MIN and MAX of [U Algorithm 1]. The subphase ends when a DDES discovers the 
first augmenting path. At that point the Augmentation Subphase begins. 

The Assignment Subphase is based on a search forest F. Section 12.11 defines F and discusses 
its properties. Then Section [2.21 gives the algorithm for the subphase. 
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2.1 The Search Forest 


At any point in a phase consider the following forest F (which represents the current search struc¬ 
ture). F has a node for every bd* v that has a finite minlevel. Each free petal* is a root. (At the 
start of the phase this means each free vertex is a root.) The parent of v corresponds to the first 
predecessor of v. More precisely suppose MIN scans an edge {u, v) and changes minlevel(u) from 
oo to a finite value. Here u is a vertex, and bud*(ri) is a node of F. The new tree F has bud*(ri) 
the parent of v. Thus the F-path from a node v (a bud*) to its root, say F{v), corresponds (by 
contraction) to a minlevel(u) path of props from u to a free petal*. 

This completes the definition of F. Now we describe how F gets modified by the MIN and 
MAX procedures. 

If MIN assigns a first predecessor to v, F gets a new leaf v. 

Suppose MAX forms a new petal P for bridge rg of tenacity t = 2i + 1. Let P* be the 
corresponding petal*, and let b be its bud (i.e., the bottleneck of the DDES for rg). Take any 
V G P* that is a bd* at the start of the current DDES. Lemma ll.lH l shows the u6-subpath of F{v) 
is contained in P*. So Lemma ll.ll iiii) shows the entire petal P* corresponds to a subtree S of F. 
We conclude the new F is formed by contracting the edges of S' in F. 

2.2 The Algorithm 

We use incremental-tree set merging as the data structure for F. Let D be the forest of this data 
structure. A bd* 6 in F corresponds to a subtree of D that has been merged into its root node b. 
Beyond that the two forests are identical. 

The MIN step performs a grow operation on D each time a first predecessor is processed. 

A MAX step will execute DDES and either find the next petal* or find the first augmenting 
path. 

First suppose DDES, say for bridge rg, identifies a petal* F* with bud b. The DDES has 
identified the support of rg (in the trees R and G of 0). The algorithm performs a union operation 
on D for each support node. This merges the entire petal P* into its bud b. (To verify this we 
remark that a petal* Q* whose bud* q is in the support of rg is contained in F*, since every vertex 
of Q* gets a new bd* b. So the unions correctly contract F*.) The updated D corresponds to the 
new search forest F by the description of MAX in Section 12.11 

It remains to describe how the DDES is implemented, and in particular how it detects the 
first augmenting path. Assume the DDES begins with the sets of D correctly corresponding to 
the current petal*s. The DDES uses find operations on D to jump from a vertex to its bd* x; it 
moves from x using a prop. Observe that the above procedure to update F and the data structure 
D occurs after the DDES. So all the find operations in the DDES give the correct bud*. Thus 
the DDES works correctly, discovering the next petal* (which might be a free petal*) or the first 
augmenting path. 

2.3 Time Bound 

The time for the Assignment Subphase is 0{m). In proof, the time for a Micali-Vazirani phase is 
0{m) plus the time for set merging. Incremental-tree set merging uses 0{n + m) time for m find 
operations, n make-root or grow operations, and n union operations. 
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3 Augmentation Subphase 


The last search level im of the phase is implemented by the Augmentation Subphase. During this 
subphase the algorithm will know bud*(u) for every vertex v in the following sense. At the start 
of the phase bud*(u) = base*(u) from phase im-i- At some point v may enter a search stack. It 
remains in the stack until it is deleted, either because of an ap or because it is known not to be 
in any ap. From this point on bud*(u) will always be computed using the data structure for the 
stack. 

The stack models this variant of DBFS; A DBFS anchored on A, for bridge rg, starts with a 
(shortest) path A from level /q to a vertex r which is on a bridge rg of tenacity Denote A as 
the sequence Ai, A 2 , ■ ■ ■, A£, where each Aj is a vertex or a petal* (a vertex A* will be a vertex of 
finite oddlevel not in any petal*). As usual the DBFS visits some set of vertices, say X, and finds 
either an ap or a new petal* which actually is X. The anchored DBFS has this defining property: 
If an ap is found then V(A) C X. If A is a petal*, say with bud b, then b = bud*(Aj) for some i, 
and V{A) n X consists of the sets Aj,j > i. 

We note that the algorithm will be valid for a pointer machine. 


3.1 Data Structure 

Array A stores the sequence V (Ai), ... ,V (A^). Each vertex v has a pointer (index) 


a[u] 


0 V(A) 
i A\i\ = V. 


A companion array gives the buds for A: 

l?[i] = bud*(Aj), i = 

A stores bud*(Aj) before the other vertices of A*. Thus for a vertex of G in V{A), bud*(u) is the 
vertex b = B[i] where i is the maximum index satisfying a[b] < a[u]. 

In addition there are some bookkeeping lists, that are easily maintained: Each petal* P has a 
list of props uv, u G P and u a predecessor of v, and bridges uv of tenacity im, u G P {v may or 
may not be in P). Each blossom (from im-i) has a list of its G-vertices (it will be used to populate 
A). 


3.2 Algorithm 
Highlevel Outline 

A is managed as a stack. The main loop computes A and then does a DBFS anchored on A. 
If an ap is discovered all of V{A) gets deleted from the graph (by definition of anchored DBFS). 
If a bottleneck b is discovered, b = bud*(Aj), we can merge petal*s Aj,..., A^ (again definition 
of anchored DBFS). So pop the B stack to erase the buds above b, and add the new vertices of 
the petal* to A. Also merge appropriate bookeeping lists for the new petal*. Then do a dfs from 
Ai, using props, trying to complete a new A-path (i.e., search for a base or blossom that has an 
incident bridge). If one is found start the next iteration. (Do the same thing in the ap case.) If 
not, Aj cannot be on any ap. Delete Aj from the graph, and dfs search from Aj_i. 

Detailed Algorithm 
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3.2.1 Computing A 

In the following, each time A is extended by a new vertex or petal* (which must be a blossom plus 
its base) Ai, push V{Ai) onto A and push bud*(Aj) onto B. 

We start with A as the sequence Ai,...,Ai {i = e) from the previous search; if A is empty, 
initialize it to a free petal* Ai, and set i = 1. 

Do a dfs from Ai, using props to find the next vertex or petal*. Extend the dfs path A (0 
or more times) until it reaches a vertex or petal* Ai that either is either a dead end or it has an 
incident bridge. 

If Ai has no outgoing props, it cannot be in any ap. Delete Ai from A and the graph. Continue 
the dfs. 

If Ai has an incident bridge rg, do a DDES anchored on A. 

3.2.2 Anchored DDES 

At any time in the DDES say each DFS path ends at its active vertex (called “center of activity” 
in IS]). We will maintain the invariant that some vertex a G V{A) is active; more precisely there 
is always an index e such that a = bud*(Ae) is active and the DDES has visited bud*(Aj) iff i > e. 
Whenever the algorithm moves a, it descends to the preceding vertex of A. Thus the invariant 
is preserved. The invariant implies that the DDES correctly implements anchored DDES (i.e., if 
an ap is found then V{A) C A; if not then A is a petal* with bud bud*(A£) for some e, and 
V{A) n A = Ui>e ^(^i)). Let u denote the other active vertex (perhaps u = a). 

The algorithm always records e. It moves a by decrementing e by 1. 

The rule for how the DDES moves is determined by its state. The 2 possible states are called 
Disjoint (when u ^ V{A)) and Meet (when u = a). (There may be > 1 move between states, e.g., 
the algorithm starts backtracking when it enters state Meet.) 

State Disjoint, u (A). 

If level(a) > level(tt), move a to the preceding vertex of A. Stay in state Disjoint and make the 
next move. 

If level(ri) > level(Q:), let uv be a prop from u to a lower level. (It exists since we’re not 
backtracking.) If a[u] = 0, move u to bud*(u). Stay in state Disjoint and make the next move. 

If not (i.e., a[u] / 0) then u G F(A). If e > 1 let 6 = B[e — 1], i.e., the bud preceding a in A. If 
a[v\ > a\b], bud*(u) has been visited, so reject prop uv. Otherwise a[u] < a[b], so a[u] < a[a\. Move 
a 0 or more times until a[Q!] < a[u]. Move u to a. Now the 2 DFS paths have met so enter state 
Meet. 

State Meet, u = a. 

Enter backtracking mode. Proceed as usual. Note that a always remains active during the 
backtracking. As before the backtracking DES always backs off from a previously visited vertex. 
It scans the next prop following the rule for state Disjoint, until (as usual) it finds a vertex u ^ a 
at level < level(Q:). If n ^ ^(^) h goes to state Disjoint. If tt G F(A), it will reenter state Meet at 
the new, lower a. 

Initializing the DBFS for bridge rg. 

Set e = I, a = bud*(Ae). If a[g\ = 0 then set u = bud*( 5 ') and begin the DDES in state Disjoint. 
Otherwise g G F(A). If a[g] > a[ot\ then r and g are in the same petal* so the DDES terminates. 
Otherwise a\g\ < a\a\, i.e., <7 is in a petal* of A preceding A^. As in state Disjoint, move a to catch 
up with g, and enter state Meet. 
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3.2.3 Finishing the DDFS 


If an ap is discovered, delete it from the graph, and (as in [2j), repeatedly delete all vertices and 
petal’s with no predecessors. This deletes V{A), by definition of anchored DDFS. 

If a bottleneck b is discovered, b = a = bud*(Te). Merge petaDs In the data 

structure pop the B stack to erase the buds above b, and add the new vertices of the petal* to A. 
Also merge appropriate bookeeping lists into those for the new petal*. 

Now in either case start the next iteration. 

3.3 Time Bound 

An anchored DDFS is a special case of DDFS. So as in the Micali-Vazirani algorithm, all the 
anchored DDFSs use total time 0{m) plus the time for set merging. It is easy to see we use 0(1) 
time to identify the bd* of a vertex in A. Thus all anchored DDFSs use total time 0{m). 

Beyond DDFS the algorithm performs a number of dfs’s to extend A. An edge of the graph 
gets scanned at most once in all these searches. Hence they use total time 0{m). 
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